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Abstract 



> 

, We perform a calculation of the isospin and hypercharge axialvector current 

O ! propagators {Aj^^{q) and A^g(g)) to two loops in SU (3) x SU (3) chiral pertur- 

bation theory. A large number of 0{p^) divergent counterterms are fixed, and 
Q^ ' complete two-loop renormalized expressions for the pion and eta masses and 

decay constants are given. The calculated isospin and hypercharge axialvec- 
tor polarization functions are used as input in new chiral sum rules, valid to 
Oh! second order in the light quark masses. Some phenomenological implications 

pH ' of these sum rules are considered. 

> 

[X^ ■ I. INTRODUCTION 

Although low energy quantum chromodynamics remains analytically intractable, the 
calculational scheme of chiral perturbation theory [|T| (ChPT) has led to many valuable 
contributions. Following the seminal papers of Gasser and Leutwyler 0,^ , numerous studies 
conducted in the following decade convincingly demonstrated the power of ChPT. The 
state of the art up to 1994 is summarized in several reviews e.g. |QJ^ (see also 10,0). The 
exploration of ChPT continues to this day, and two-loop studies represent an active frontier 
area of research. These include processes which have leading contributions in the chiral 
expansion at order HHUI even , as well as systems for which precision tests will 



soon be available, e.g. the low-energy behaviour of vrvr scattering |T^-|T^. While the case of 
SU (2) X SU{2) ChPT to two-loop order has been relatively well explored (in particular see 
0), works in SU{3) x SU{3) ChPT are still few in number [0,|rT[|T|,0|ril . 

Recently, we performed a calculation of the isospin and hypercharge vector current prop- 
agators (AY3(g) and AYg(g)) to two-loop order in SU{3) x SU{3) chiral perturbation the- 
ory 0. A partial motivation for working in the three-flavour sector stems from its inherently 
richer phenomenology. In particular, it becomes possible to derive new chiral sum rules which 
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explicitly probe the 5'f/(3)-breaking sector. With the aid of improved experimental informa- 
tion on spectral functions with strangeness content, it should become possible to evaluate 
and test these sum rules. 

We have completed this program of calculation by determining the corresponding isospin 
and hypercharge axialvector current propagators (A^3(g) and A^g(g)). Determination of 
axialvector propagators is much more technically demanding than for the vector propagators, 
but at the same time yields an extended set of results, among which are: 

1. a large number of constraints on the set of 0{p^) counterterms, 

2. predictions for threshold behaviour of the Svr, KKtt, KKtt, rjinT, etc axialvector spec- 
tral functions, 

3. an extensive analysis of the so-called 'sunset' diagrams, 

4. new axialvector spectral function sum rules, 

5. a new contribution to the Das-Mathur-Okubo sum rule [^, and 



6. a complete two- loop renormalization of the masses and decay constants of the pion and 
eta mesons. This final item places the axialvector problem at the heart of two-loop 
studies in SU{3) x SU{3) ChPT. 

We begin the presentation in Sect. II by presenting basic definitions and describing the 
calculational approach. To illustrate the procedure, we summarize the results for the tree- 
level and one-loop sectors. Our calculation of the two-loop amplitudes and the corresponding 
enumeration of the 0{p^) counterterms is given in Sect. III. The construction of a proper 
renormalization procedure forms the subject of Sect. IV. It provides the framework for the 
removal of divergences, described in Sect. V, and leads to finite renormalized expressions for 
the meson masses, decay constants, and polarization functions. We give explicit expressions 
for the isospin polarization functions in Sect. VI and continue the presentation of finite 
results in Appendix B. Sect. VII deals with the determination of spectral functions. The 
subject of chiral sum rules is discussed in Sect. VIII and Sect. IX contains our conclusions. 
Various technical details regarding sunset integrals are presented in Appendix A and as 
mentioned, our final expressions for meson masses, decay constants and the hypercharge 
polarization amplitudes are collected in Appendix B. In addition, at several points in the 
paper we compare results as expressed in the 'A-subtraction' renormalization used here with 
a variant of the MS scheme. 



II. TREE-LEVEL AND ONE-LOOP ANALYSES 

The one-loop chiral analysis of the isospin axialvector-current propagator was first car- 
ried out by Gasser and Leutwyler who used the background-field formalism and worked in 
an SU (2) basis of fields 0. We shall devote this section to a re-calculation of the isospin ax- 
ialvector propagator through one-loop order, but now done within the context of a Feynman 
diagram calculation and using an SU{3) basis of fields. 
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A. Basic Definitions and Calculational Procedure 



Our normalization for the SU (3) octet of axialvector currents is standard, 

^^ = 5"y7^5g (fc = l,...,8) . (1) 
In this paper, we shall deal with the axialvector current propagators 

^ZiQ) = ijd*x e"^-"" (0|T (^^(a;)A^(0)) |0) (a = 3, 8 not summed) , (2) 
having the spectral content 

-Im AZiq) = {q'q'' - q'gn&iq") + , (3) 

TT 

where and arc the spin-one and spin-zero spectral functions. This motivates the 
following tensorial decomposition usually adopted in the literature, 

Kaiq) = W - q'rn^^aiq') + q'^q^'^faiq') , (4) 

where 11^^^ and 11^^^ arc the spin-one and spin-zero polarization functions. The low-energy 
behaviour of the spin-zero spectral function p^l is dominated by the pole contribution as- 
sociated with propagation of a Goldstone mode, 

pflis)^ Fl 5{s - Ml) + p2{s) . (5) 

In the chiral hmit, one has — > and PaI{s) — > as well. 
The lowest-order chiral lagrangian 

£(2) 

is given by 

£(2) ^El^ (d.UD^U^) + ^ IV {xU^ + Ux^) , (6) 

where Fq is the pscudoscalar meson decay constant to lowest order and x = ^BqIU is 
proportional to the quark mass matrix m = diagonal (m, m, m^) with 

Bo^-^{qq) . (7) 

Note that we work in the isospin symmetric limit of m„ = = m and that to lowest order 
we may use the Cell Mann-Okubo relation, 

ml = l{^ml-ml) . (8) 

The field variable U is defined in terms of the octet of pseudoscalar meson fields {0jk}, 

U = exp{iXk ■ (l)k/Fo) , (9) 
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and we construct the covariant derivative D^U via external axialvector sources 

D^U = d^U + ia^U + iUa^ . (10) 
The axialvector source has a component for each of the SU (3) flavours, 

= < • (11) 

Adopting the approach carried out in Ref . , we make use of external axialvector sources 
to determine the axialvector propagators. The procedure is simply to compute the S- 
matrix element connecting initial and final states of an axialvector source. Analogous to the 
analysis of Ref. [§| , the invariant amplitude is then guaranteed to be the axialvector-current 
propagator, say of flavour 'a', 

{a,{q\\')\S-lKiq,\))=^i2nr5^'\q'-q) e^q' , X')AZ{q)eA<l, ^) ■ (12) 

We shall typically use the invariant amplitude symbol J^^^a to denote various individual 
contributions (tree, tadpole, counterterm, IPI, etc) to the full propagator. 



B. Tree- level Analysis 

For definiteness, the analysis in the remainder of this section will refer to the isospin 
flavour. Given the lagrangian of Eq. (^, it is straightforward to determine the lowest-order 
propagator contributions, 

= F^g,. - , ^i^. q,q. . (13) 

The two terms represent respectively contributions from a contact interaction (Fig. 1(a)) 
and a pion-pole term (Fig. 1(b)). Although M.^^^^^ has an exceedingly simple form, it is 
nonetheless worthwhile to briefly point out two of its features. First, as follows from unitarity 
there is an imaginary part corresponding to the pion single-particle intermediate state, 

TmM%t =7^F',6{q'-ml)q,q, ■ (14) 



However, there is also a non-pole contribution to M.^i^l'^^^ ■ Its presence is needed to ensure 

the proper behaviour in the chiral limit {df^A^ = 0), where TZe Ai^j^H^^ is required to obtain a 
purely spin-one (or 'transverse') form. This is indeed the case, as we find by taking — >■ 
in Eq. dg). 



^2 

{q,i(iu - q^g^u) ■ (15) 



Throughout this paper we adopt the phase employed in Ref. |21|, which is opposite to that used 
in Ref. [i. 



C. One-loop Analysis 



At one-loop level, the axialvector-current propagators are determined from the lagrangian 
of Eq. and correspond to the Feynman diagrams appearing in Fig. 2. The loop correc- 
tion to the lowest-order contact amplitude appears in Fig. 2(a) and Figs. 2(b),(c) depict 
corrections (whose significance we shall consider shortly) to the pion-pole amplitude. We 
find for the one-loop ('tadpole') contributions to the isospin propagator, 



_^(tadpoie) ^ _^ ^2A(m:) + A{ml)) g,, + 



Ai 2A{ml) + A{m\) 
3 — 



+ 



i A{ml){ml - Aq^) + 2A{m]^){m\ - q^) + A{ml)ml 
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(g2 



2^2 



(16) 



where all masses occurring in pole denominators are understood to have infinitesimal nega- 
tive imaginary parts. In the above, A is the scalar integral 



A(m^ 



dk 



/j.2 _ ^2 



(17) 



where dk = d'^k/{27iY. Hereafter, any integration measure accompanied by a super-tilde 
will have a similar meaning. The evaluation of A is standard, and we have 



A{m^) 



fx 



A-d 



A-d 



r 1 



d-4 



-2im^X 



1— im 



[m 



2\(d-2)/2 



167r' 



with 



d-4 



167r2 



^-^(log4vr-7 + l) 



(18) 
(19) 

(20) 



The quantity introduced in Eq. (|T8D is the mass scale which enters the calculation via the 
use of dimensional regularization. The /^"^"^ prefactor in Eq. (p!9|) ensures that A[m?) has 
the proper units in d-dimensions. 

To deal with divergences arising from the loop corrections, one must include counterterm 
amplitudes. It suffices to employ the well-known list of counterterms {Li] {i = 1, . . . , 10) 
and {Hj} (j = 1,2) appearing in the 0{p'^) chiral lagrangian of Gasser and Leutwyler 
Analogous to Eq. (|I9|) for the yl(m^) integral, each 0{p'^) counterterm is expressible as an 
expansion in A, 



(d-4) 



(0) 



(21) 



n=l 



^See also the discussion surrounding Eq. (26) of Ref. 
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where the leading degree of singularity is seen to be linear. The counterterm diagrams involve 
a contact term (Fig. 2(d)) as well as contributions to the pion pole term (Figs. 2(e),(f)) and 
yield the following isospin counterterm amplitude, 



M 



(CT) 



2(Lio - 2Hi){q^q^ - q^g^^) + 8 {ml + 2ml) + ^IL. 



(ml + 2m%)LA + miLr. 



q- 



mt 



^ q^ {{ml + 2m\)U + m^^s) - 2ml (« + 2m^)^6 + ^^^s) 

[q^ — m^)^ 



(22) 



Results through One-loop Order 

Our complete expression for the isospin axial vector current propagator through one- loop 
order is given by the sum of Eqs. 



A 



, .(tadpole) A^(CT) 



(23) 



The resulting expression is complicated and seems to lack immediate physical interpretation 
because we have not yet accounted for renormalizations of the pion's mass and decay con- 
stant. A detailed account of the renormalization procedure is deferred to Sect. IV. However, 
we note here that the renormalized masses and decay constants have the expansions 



p2 ^ ^2(0) _^ p2(2) p2{4) _^ _ _ _ 



(24) 



where we have temporarily suppressed flavour notation and the superscript indices {(«)} 
denote quantities evaluated at chiral order {p*}. To one loop 0, the explicit expressions are 
given hyf\ 



Mi 



f2 
-^0 



ml + 



[ml + 2ml)Lf + m^L, 



2 r(0) 



2ml_ ml 



12- 



167r2 



12- 



4 2 

mi , mt 
^In^ 



2 2 2 

m^m^ m^ 
■ In — 7- 



327r2 /i2 96n- 

(0)^ , _2/r(0) or(0)- 



-8ml {{ml + 2ml){Lf - 2Lf) + ml{Lf^ - 2L[ 



(25) 



Upon combining the information gathered in Eqs. (p3D,(p5D as well as the divergent part of 
the one-loop functional given in Ref. f^, one finds through one-loop for the renormalized 
isospin propagator. 



^A3 



Flg^^^ + 2(lS? - 2H^^^){q^q^ - q^g^^^ 



Ml 



(26) 



one-loop order, our counterterms {Lf*^} and {H^^^} are equivalent to the {L^} and {H^} of 
Ref. il. 
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The hypercharge propagator can be obtained analogously. We comment that both 
the isospin and hypercharge amplitudes contain the regularization dependent constant h[^\ 
and are therefore unphysical. A physically observable quantity is obtained from the differ- 
ence, 



^A3 



^A8 



9'" - q'q" 



V 



q2 _Ml g2 _ 



(27) 



III. TWO-LOOP ANALYSIS 

The general structure of propagator corrections is displayed in Fig. 3, which includes the 
one-particle irreducible (IPI) diagrams of Fig. 3(a) and the one-particle reducible (IPR) di- 
agrams of Fig. 3(b). The latter consists of both vertex and self-energy corrections. Through- 
out this section, we continue to focus on the isospin sector when giving exphcit expressions 
for two-loop amplitudes. 



A. Two-loop Analysis: IPI Graphs 

First, we consider the IPI diagrams of Fig. 4 and Fig. 5. The graphs of Figs. 4(a)-(c) are 
simple in the sense that they are either equal to or the negative of identical graphs in which 
external vector sources occur. As such, they can be read off from the work of Ref. 0. This 
is only partly true of Fig. 4(d), and the graph of Fig. 5 (the so-called 'sunset graph') has 
no counterpart in the vector system. For convenience, we compile definitions and explicit 
representations for the sunset-related functions in App. A, leaving detailed derivation of 
these results for another setting [p^] . 

The amplitudes of Figs. 4(a)-(b) are of the general form 

M,A^aAh] = 1^ E <^tt^m , (28) 

where the {a^'^''''*} are numerical coefficients (in some cases dependent on dimension d) and 
the {A{k)} are the integrals of Eq. (|17|). The sum over the indices k,C. simply reflects the 
need to perform flavour sums independently for each of the two loops. The amplitude of 
Fig. 4(c) can be expressed similarly, 

M,A^c] = E U , (29) 

except now the coefficientsja^'^^} are proportional to squared meson masses and there are 
{Li} factors arising from the 0(p^) lagrangian. 

Each of the above amplitudes will diverge for d — > 4, and we write 

M,.s[4a, 46, 4c] = |^ la^^^'^^^X' + a^^^'^^^X' + . . .1 , (30) 
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{a,b,c) 



are 



where the singular quantity A has been previously defined in Eq. (^Oj), 02"'"''''' and 

numerical quantities, and the ellipses refer to terms which are nonsingular at c? = 4. As 

2 

expected for two-loop amplitudes, the leading singularity goes as A . 

For Fig. 4(d), the contribution from Liq is the negative of the vector case, but there are 
also many new terms, 

4i 



-Mm-sN = -{q^q, - g^,q')-;^i2A{7r) + A{K))L 



10 



^9t^u 



32 



64, 



A(7r) [mi[{A8 + — )Li + (16 + — )L2 + (24 
- 8(5m2 + 4m^)L4 - 28m^Lf 



64 16 
+A{K) [ ml[QAL, + —L2 + (16 + — )L3 

- 8iml + Qm\)U - 8{ml + 



(31) 



K 



16 



8 16, 



+A(r/) m^[16Li + —L^ + (- + —)L^] - 8m' U - i^m^L 



d 



3d' 



Aside from the presence of a transverse component proportional to {q^Mqu—QpLvq^), it resembles 
the form of Fig. 4(c) and shares the degree of singularity shown in Eq. (pO]). 
Finally, there is the sunset contribution of Fig. 5, which we write as 



= ^^fiuiq'^, ml, ml) + ^Hf^uiq^, m\, m\) 
+ T^^M^'(9^ "^l^ ^k) + ^f^uiq^, ml, m\) 



(32) 



The quantities Ti^jy and £^1, are defined in App. A by Eqs. ( |A1| ),([A^. The fiavour depen- 
dence of the individual sunset contributions can be read off from the arguments of the above 
functions. The defining representations for 7Y^,y and (c/Eqs. ([TO]), ([7^)) are given as 
multidimensional Feynman integrals. It takes considerable analysis to reduce these inte- 
grals to forms which can be compared directly with the other two-loop amplitudes, and one 
finds that the sunset amplitudes share the singular behaviour of Eq. (pO). 



B. Vertex Corrections 



The IPR graphs are of two types, the vertex modifications of Figs. 6(a)-(e) and the 
self-energy effects of Figs. 7(a)-(e). 

The two-loop vertex amplitude F^^^^ is defined by 

{P,{q')\S - iMq, A))Fig. 6 = 7(27r)V^)(g' - q) e^{q, A)rg) 
A typical example is the vertex amplitude of Fig. 6(b), 



(a = 3, 8) . 



18F3 



A{ti) 16A(7r) + 12m^5(0, tt) + llA{K)A{'K 



(33) 
(34) 



^QA^{K) + A{r]) ?,A{K) - AmlB{{], tt) + Am\B{{], K) 



where _B(0,m^) is defined and evaluated as 



5(0, 



dk 



A{im?) - - — -A{m^) 



The amphtudes of Figs. 6(a),(b) have the generic form 



(35) 



(36) 



with numerical coefficients f^."'^'' whereas those of Figs. 6(c), 6(d) are 

-^0 kl 



(37) 



with coefficients f^^''^^ proportional to the squared meson masses. As in Eq. (|3^)), all the 
invariant parts of non-sunset vertex functions are found to be independent of the propagator 
momentum g^. 

For the sunset contribution of Fig. 6(e), we write 



-h^ [q ; m^; m^; m\ + —h^ [q ; m^; m^; 2(m^ + 



36 



, ^ T f 2 2 2^/2 2 \\ , ^ T f 2 2 



m 



K 



(38) 



The vector- valued quantities Ji^ and are defined in Eqs. ( |A6[ ) , (|A7D of App. A, and are 
found to share the singular behaviour of Eq. (BO). 



C. Self-energy Corrections 

The two-loop self-energy T,^^ arises from the meson-to- meson matrix element, 
{Pa{q')\S - l|Pa(g))Fig. 7 = z(27r)^5W(g' - q) Si^) (a = 3,8) . 



(39) 



The non-sunset contributions of Figs. 7(a)-7(b) are proportional to {A{k)A{i)} and 
those of Figs. 7(c)-7(d) are proportional to {A{k)L£}. As an example, the isospin self- 
energy contribution of Fig. 7(a) is 



"9" 



+ ( — ] m 

'24 9 



A'in) + 



~36 



A{n)A{r^) 



+ 
+ 
+ 



7g^ 7ml + 2m^ 2ml + "^1: 



1^ 



9 
.2 



17m|^ ^ 7{ml + 2m\) 



18 
2g2 

Ts" 45 ■ 90 

g2 Qm\ + ml ml + m'j 

30 ~ 



A{7t)A{K) 



A\K) 



(40) 



30 



+ 



45 



A{K)A{v) + ■-^A'ir^) 
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As in Eq. (|40|), the remaining non-sunset self-energies are at most linear in the propagator 
momentum g^. 

Finally, the sunset contribution of Fig. 7(e) is given by 

Fo'^? [e] = "^Siq^; ml; ml) + ^^(g^; ml; m^) 

+ ^^(^^; ^l'^ ^l'^ ^l) + 2*^"^^ + ^"k)) 

12 1 

j^_R^q^-rnl;m\;-{ml-m\)) + -U{q^;ml;m\) , (41) 

where the quantities R and U are yet new sunset functions defined by Eqs. (|A2|) , (|A8|) , (|A9|) 
of App. A. Both they and the non-sunset self-energy contributions share the degree of 
singularity of the IPI and vertex amplitudes (c/Eq. (^OD). 



D. 0{p^) Counterterms 



To construct the counterterm amplitudes needed to subtract off divergences and scale 
dependence contained in the two-loop graphs, we refer to the lagrangian of 0{p^) coun- 
terterms of Fearing and Scherer |Q. From their compilation, we extract 23 counterterm 
operators which contribute to the axialvector propagators. 

The 0{p^) counterterm amplitudes are computed in like manner to the two-loop contri- 
butions considered thus far in this section. For example, the isospin IPI two-loop amplitudes 
will require a corresponding counterterm contribution A^^,^3 as computed from 

(a3(g', A')|5(^*) - l|a3(g. A)) = i{27,Y5^'\q' - q) e^\q\ A') e^{q, A) , 

and one obtains 



(42) 



= (q^qu - q^gt^u) 



4{ml + 2ml){B29-2Big) 



Aq^iBss-2B32) + 8mliB2s-B, 



46 J 



(43) 



4{ml + 2mj^yB2i - 4(3mt - 4m>^ + 4m|^)5 



19 



Amliml + 2m^)(5i6 + B^s) + 4mt(25i4 - B^, - B 



17) 



This has three independent 0{p^) counterterms for the 'transverse' amplitude proportional 
to q^qu — q^g^u and five for the 'longitudinal' amphtude proportional to g^^. 

Employing a relation analogous to Eq. ( p3D for the vertex counterterm amplitude, one 
finds for the isospin case. 



,{ct) 
fj.3 



AmliBu - Bi-j) - 4 (4m^ - Am^ml + 3m^) B 
I + 2m\) {2ml{Bie + 2Bi^) + A{ml + 2ml)B2i 



(44) 
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and a relation analogous to Eq. (^) leads to the isospin counterterm self-energy 



^3 - J2 



ml {Wi + 2B2) + (5mt + 4m^)S3 



+ {ml + 2mj^) [2mlBi + 3m^(m^ + 2mi.)5f 



+ {ml + 2m^) (m^^is + {ml + 2m^)52i 
Finally, we express the {-B^} in dimensional regularization as0 



(45) 



n=2 



(2)/,, A t2 



>(0), 



This representation, together with Eq. ([T9|) for the A-integral and Eq. (|2T|) for the 0{p^) 
counterterms {L^}, expresses the axial vector propagator as an expansion in the singular 
quantity A. 

IV. RENORMALIZATION PROCEDURE 



The axialvector propagator A^^ will have contributions from both the IP/ part M.'^ 



and the IPR pole term. 



-ml + T.a{q^) 



(a = 3, 8) . 



(47) 



In the above, the 0{p^) counterterms are understood to be already included in M.'^'^ ■, in the 
self energy and also in Ta{q^)- The latter is defined in terms of the vertex amplitude 
r^^(g) as 



(48) 



The renormalized mass Ma and decay constant Fa are defined as parameters occurring in 
the meson pole term 



q^ -ml + T.a{q^) 



p2 

+Ra{q') 



(49) 



where Ra{q^) is a remainder term having no poles. 



^The dependence of L^"^ (/u) and B^^^ (/i) upon the scale // is determined from the renormalization 
group equations and has been explicitly given in Ref. 0. 
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A. Identification of the Meson Mass 



From Eqs. (P),®, it follows that is a solution of the implicit relation 

M! = ml - J:,{M!) . (50) 

Since we have already calculated S(m^) (in the following, we temporarily omit flavour in- 
dices), it makes sense to expand the self-energy S(M^) as 



S(M2) =S(m2) + S'(m2)(M2-m^) + ... . 
Then, expressing the squared-mass perturbatively, 

and similarly for the self-energy, we obtain the perturbative chain 



M(2)2 = m 



where we have noted that 



-SW(m2) + ... . 



(51) 



(52) 



(53) 
(54) 
(55) 



(56) 



We briefly exhibit this procedure at one-loop level. Using the fourth-order isospin self- 
energy 

i2 A^2 ^2 ^2 ^2 



m: 



K + 2mj,)(g^L4 - 2miL^) + m^(g^L5 - 2m^L; 



(57) 



we obtain from Eq. 



3A(7r) - A(r/) 
6 



5m: 



+ 2m2,)(L4 - 2U) + m2(L5 - 2U) 



(58) 



Expanding the A-integrals in a Laurent series around d = 4, we readily obtain the result 
cited earlier in Eq. (p5[) . In like manner, the one-loop hypercharge self-energy 



8 

16, 



(m 



K 



(59) 



yields 



12 



16 

y 



mliml + 2ml) - 2^6) + mjLg 



(60) 



These agree, of course, in the SU{3) hmit of equal quark mass. The analysis at two-loop 
level proceeds analogously. 



B. Identification of the Meson Decay Constant 



With the aid of Eq. (|50|), it is straightforward to write the meson pole term in the form 
(again temporarily suspending flavour indices) 



where S(g^) is the divided difference 



tiq 



q^-M^ 



From the definition in Eq. 



q^-M'^ 

of the squared decay constant F^, we have 
T\q') (l-S(g2) + (S(g2))2 + ... 



= lim 



Let us analyze this relation perturbatively. 

We begin by writng the vertex quantity r(g^) (evaluated at = M^) as 

r(M2) = r(o) + r(2) + r(^)(M2) + ... , 

where both r(°) and T^^) are independent of g^. Expanding r(^)(M^) as 
rW(M2) = r(^)(m2) + r(^)'(m2)(M2 - m^) + . . . , 



(61) 



(62) 



(63) 



(64) 



(65) 



we see from chiral counting that to the order at which we are working, one is justified in 
replacing T^^^M^) by T^^\m^). As for the self energy dependence in Eq. (|63D , we first 
observe that 



S(M2)= hm S(g2) = S'(M2) . (66) 
Recalling from Eqs. (p^) , (^9]) that T^^^^q"^) is linear in g^, we have the perturbative expression 
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(67) 



where the error in replacing ^(^^'(M^) by S(6)'(m2) appears in higher order. Thus, the 
perturbative content of Eq. (|63|) reduces to 



x(l-S(^)'-S(6)'(m2) + (sW')2)+... . 



(68) 



Upon organizing terms in ascending chiral powers, we obtain the following expression, valid 
through two-loop order, for the decay constant 



p(2) _ ^^(0)^(4)' 

2 



+ 



r(4)(^2)_ ^^(2)5.(4)' 



+ rw(-ls(^)'(m2) + ^(s(^)')2 

2 o 



(69) 



where we have collected together terms of a given order. 

As an example, Eq. ( |69D readily provides a determination of the isospin and hypercharge 
decay constants through one-loop order. The corresponding vertex quantities are 



A [2A(vr) + A{K)] + A 



[ml + 2m]^)Li + m^Lr 



+ 



2i 



r. = F,- -A{K) + - 

-To -To 

.(4)' 



2 ' 2m^)L4 + mjL5 



+ . . . 



Upon inferring {k = 3, 8) from Eqs. (p^) , (pUD , we obtain 



(70) 



i.(0) + ^(2) = Fo - — [2A(vr) + 

[ml + 2m^)L4 + m^Lg 



4 



2Fo Fq 



[ml + 2m|')L4 + m^Lg 



(71) 



The one- loop expressions for the pion decay constant (appearing in Eq. (pSD) and the eta 
decay constant (not shown) follow from the above relation. The two-loop corrections are 
found analogously. 



C. The Remainder Term 



The preceding work allows extraction of the remainder term R{q^) defined earlier in 
Eq. (^). Making use of Eqs. ([5^)-(|55D as well as Eq. (pUj), a straightforward calculation 
yields the expression 



2Ff 



q2 _ ^2 



^2 S(6)(g2) _ S(6)(m2) - (g2 _ m2)S(6)'(m2) 



(g2 — m'^y 



(72) 
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It is manifest that R{q'^) has no pole at = m^. Moreover, since the non-sunset vertex 
functions F^^^^ [a]-[(i] of Sect. 3 are constant in g^, they do not contribute to R{q^). Nor do 
the non-sunset self-energies S*^^) [a]-[(i] since they are at most linear in q"^. Thus, the only 
contributors to R{q^) are sunset amplitudes, and it is straightforward to obtain the isospin 
and hypercharge remainder functions directly from Eq. ([7^). 

D. The Polarization Amplitudes 11^^^ and n^'^) 

From the tree-level and one-loop results, we anticipate that the two-loop IP/ amplitude 
Ai^ff} will contain an additive term involving the meson decay constant, and can thus be 
written 

Ml'^^g.AF'P+M,. , (73) 

where Xi^u denotes the residual part of the IPI amplitude. It turns out that much of the 
rather complicated content in the IPI amplitudes of Sect. 3 is attributable to the two- loop 
squared decay constant (F^)'^^\ as can be verified from decay constant results derived earlier 
in this section. First, we re-express the expansion of Eq. (|68D as 

+ + + , (74) 

with (F2)(o) = {F^Y^^ = 2FoF(2) and 

(F2)W = (F(2))2 + 2FoFW . (75) 

One then compares the gfj,^ part of with {F'^Y'^\ The residual amplitude A4^„ is simply 
the difference of these. 

It is convenient to replace the residual amplitude Ai'^'^ and the remainder term 
—q^q'^R{q^) (which arises from the meson pole term but itself contains no poles) by equiva- 
lent quantities n^^-* and 

>l^'^(g)-g'^g'^i?(g2) = (gMg--gV')n^'^(g')+^7'''n(°)(g2) . (76) 
We shall employ Fl*^^^ and Fl^'') throughout the rest of the paper. 



V. REMOVAL OF DIVERGENCES 

Thus far, we derived lengthy expressions for the various two-loop components of the 
axialvector propagators and then determined the renormalization structure of the masses 
and decay constants. At this stage of the calculation, there are many terms which diverge as 
(i — i> 4 and which must therefore be removed from the description. Below, we carry out the 
subtraction procedure by expanding the relevant quantities in powers of the parameter A and 
then using 0{p^) and 0{p^) counterterms to cancel the singular contributions. In particular, 
this process will determine a subset of the so-called /3-functions of the complete 0{p^) 
lagrangian. This is of special interest because the divergence structure of the generating 



functional to two- loop level can be obtained in closed form and our results derived in 



the following can be used as checks of such future calculations. 
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A. Removal of Dependence 



It will simplify the following discussion to define the counterterm combinations 

A = 2Bu — Bn , B = Biq + Bi^ , C = Bi^ — B20 , 

D^BiQ + B2i, E = BiQ -B21, F = 3Bi + 2B2 . (77) 

First we consider the decay constant sector. Upon demanding that all dependence 

vanish, we obtain six equations containing five variables. There are six equations because 
the pion and eta constants each have explicit dependence on m^, m^m\ and factors 
and the singular behaviour must be subtracted for each of these. 

We find the equation set to be degenerate, and one obtains just the following four con- 
ditions, 

^(2) _ 3^(2) ^ _ , 5(2) _ 2^(2) ^ ^3 



(J{2) ^ ^(2) _ ^(2) ^ /ygN 



24:Fi ' 72F^ ' 

13_ ^(,) ^ _73_ 

ISFo^ ' U4F^ 

The information contained in the above set is unique and any other way of expressing the 
solution must be equivalent. 

In the mass sector, we find subtraction at the A^ level to yield seven equations in eleven 
variables. The number of equations follows from the dependence of each mass on m^, m^m\ 
and m^mj^ factors (this implies six equations), along with the fact that dependence is 
absent from M^^"^. This latter fact arises becauses there is no to|- counterterm contribution 
in M^^)2, and thus there must be a cancellation between sunset and nonsunset numerical 
terms. Such a cancellation occurs and constitutes an important check on our determination 
of the sunset contribution. The equation set for the M^^^^ masses is found to be degenerate 
and just five constraints can be obtained, e.g. 



B 



f = ^ - iF(2) - 5f ) - 5(2) _ 3^(2) 
3 27Fn2 6 457 





= + + ) + + B? 

^ 81Fi 18 3^3^^ 



-B14 



(2) 

B 



— -2 + 72Sf ^ + 72Sf ^ + 216Sj^^ 



1 

48 

(2) _ JOl^ 1^(2) _ ^(2) _ 2^(2) _ 3^(2) 



216^0^ 3 

= + ^F(^) + 25f ) + ) + 3B?^ . (79) 

Finally, removal of A^-dependence for the polarization functions fl*^'^-' and fl^^^ yields a 
final set of constraints at this order, 

p(2) _ p(2) _ . p(2) _ ^„(2) 

p(2) r(2) _ 3 (2) 00(2) _ 1 

B28 - -t»46 - ' -^29 - ^^49 - ■ 
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B. Removal of A Dependence 

In a similar manner, we obtain constraints for the {B^^^} counterterms. We list these 
below without further comment, beginning with those following from decay constants. 



^(1) _ 3£;(i) 



175 



28 



(0) 



92167r2 



25 



(0) 



, 26 (0) 
+ 3 ^ 



(0) 



124°) + 124°^ 



19 



15367r2 



32^(0) 
9 ^ 



9 ^ 



-L 



(0) 



+ 



9 4 g 5 6 



Fi 



691 



829447r2 9 



(0) 



9 



Sir 



9 

44°^ - 64 



+ 



59 
18' 



(0) 



43 



30727r2 9 ^ 9 ^ 18 ^ 



34.(0) .(0) 



4L, 



(0) 



2L 



(0) 



(81) 



then masses. 



(1) 



d(1) 



2167r2 



, 40 (0) 
+ y^6 



+ 404^^ + ^4 



- -F(^) - Bi'^ - B« - 34^) , 
6 



648 



Fi 



1 



_ 4 _ 31684°) - 244°) 



63364°^ - 9792L^' 



(0) 



32164°M + 36F(^) + 216Sf ) + 216B^^^ + 6A8B^^^ 



167 



46087r2 



+ 84°^ 



644°) - ^l: 



62 



(0) 



(1) 



-B 



(1) 



371 



_ + 24°) - 1^4°) + 724°) + 244°) 

207367r2 ^ 3 ^ ^ ^ 



_ - 5f ) - 25« - 35(^) , 

3 



9 



(0) 



152 



l(°) - -L 



62 



5127r2 



(0) 



+ ^F(^) + 2Sf ) + S« + 3S« , 



(82) 



and finally from the polarizations n*^°) and 



17 



o(l) 



49 1 
576 16^^ ' 



B 



(1) 

13 



173 1 

5184 Wtt^F^ ' 



2B 



(1) 

32 



(1) 

28 



B, 



(1) 

46 



64 167r^ 



B. 



(1) 

29 



2B 



(1) 

49 



5. 



(1) 
33 



64167r2Fo2 ' 



17 1 
96 167r2 



L 



This completes the subtraction part of the calculation. 



C. A-Subtraction and MS Renormalization Schemes 



(83) 



The renormalization procedure employed originally in Ref. ||^ and adopted in this pa- 
per amounts to A-subtraction, c/Eqs. (pT]) , (|46|) . Alternatively, one could employ minimal 
subtraction (MS), 



Lt{d) 
Beid) 



/i 



2uj 



(47r)2 

,4:U> 



2uj 



(47r) 



B 



(2)MS 



{2uj) 



+ 



B 



(1)MS 



2uj 



(0)MS 



(/i,u;) + . . . 



(84) 



where uo = d/2 — 2, or modified minimal subtraction (MS), 



U{d) 

Bi{d) 



(47r)2 

(^ 
(47r)4 



2uj 



+ L^t^(/i,c.) + ... 



5 



(2)MS 



(2^) 



+ 



B 



(1)MS 



+ (4vr)X)^^(/.) + ... 



2u 



(85) 



where we make the standard ChPT choice 



1 



lnc= -^[l-7E + ln(47r)] = -C 



(86) 



Of course, there must be only finite differences between these three procedures, amounting 
to additional finite renormalizations. 

As an illustration, we relate the {B^^^} and {i^^"^} renormalization constants employed 
here to those defined in the MS approach of Ref. fl^. In the latter scheme, one writes 
further that 



LfH^^, CO) = Lr,«(/i, 0) + (/x, 0) ^ + . . . 



'MS 



MS 



and sets 



(87) 



(88) 



The 'convention' established by Eq. (^) is allowed because it can be shown [T^ that the 
effect of the quantity {fi, 0) is to add a local contribution at order p^, which can always 
be abosrbed into the couplings of the 0{p^) lagrangian. Comparison of the two approaches 
yields 
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and 



p(2) _ p(2)MS 



p(l) _ ^ p(l)MS /Q^x 

We stress that the content of Eqs. ([89|),(p0|) is partly a reflection of the convention of Eq. (p8[). 
Finally, one can combine the relations of Eq. (^) to write 

B?M = srM+(5^^Bj"-(glflf' . (91) 

We shall return to the comparison between the A-subtraction and MS renormalizations 
at the ends of Sect. VI and of App. B. 



VI. THE RENORMALIZED ISOSPIN POLARIZATION FUNCTIONS 

Having performed the removal of and A^ singular dependence from the theory, we are 
left with the A° sector in which all quantities are finite. We shall express such contributions 
entirely in terms of physical quantities by replacing the tree-level parameters m^, m|-, m^, Fq 
with M^, M|-, M^, Ft,. Any error thereby induced would appear in still higher orders. 

For any observable O {e.g. O = flli \ M^, F^,, etc) evaluated at two-loop level, there will 
be generally three kinds of finite contributions, 

O = + 0CT + OyZ . (92) 

Oj-era rcfcrs to the finite A'' 'remnants' in {A{k)A{i)} or {A{k)L£} contributions and arise 
from either the product of two A'^ factors or the product of A^ and A~^ factors. Oqt denotes 
any term containing the {Bf'^} p^-counterterms, whereas Oyz represents contributions from 
the finite Y, Z integrals (c/ Eqs. ( |A20| ),( [A2T| ) of App. A) which occur solely in sunset 
amplitudes. 

Before proceeding, we address a technical issue related to the presence of L^f^^ terms 
appearing in the 'remnant' amplitudes. Such terms are always multiplied by polynomials in 
the quark masses and hence can be absorbed by the 0{p^) counterterms, as expected from 
general renormalization theorems . In the vector current analysis of Ref . ||^ , we defined 



the following dimensionless quantities (now expressed in terms of the {Bg}), 
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i?^2F^^i?(?-(4-n4o'^) . (93) 

and thus removed all explicit l\ dependence. We repeat that procedure here by defining 
axialvector quantities 

Qa = 2F^ (--828^ + -046^) + , 

R^^F^{-Bi^ + 2B^^)+L[,'^ , (94) 

such that the counterterm dependence of n^]g (c/Eqs. (P7D,(|B^)) is identical to that estab- 
lished originally for nv3,8 ID- 

Since our finite results are quite lengthy, we restrict the discussion in this section to just 
the isospin polarization functions. Expressions for all the other observables (masses, decay 
constants and hypercharge polarization functions) are compiled in Appendix B. 



A. The Spin-one Isospin Polarization Amplitude 

Turning now to the isospin transverse polarization amplitude, we have for the remnant 



piece. 



+ log^'- 
/i 



Ml 



49 



C 



13824 



+ log 



Mi 



^288 
^ 1 



r K'J) 



192 

1 . Ml 

2~ 



1 



Ml 



/i2 \^576 1536 
5 17 C 



log 



+ - 



36864 
283 



+ 



3072 
3 C 



+ log 



768 

f^' J 

Ml 



log 



294912 4096 3072 



log 



17 
3072 

Ml 
^2 



1 , Ml' 



768 



^ r(0) 



12288 



log 



1024 

Ml 



log 



Ml- 



(95) 



where the constant C has the same meaning as in Ref. 



7e + ln(47r)] 



(96) 



Next is the counterterm contribution, 



naT(g') 



i?2 



Pa 



SMI AMI , . ^ 



(97) 



and finally upon defining 
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Wi = S-6S + 9Si , 



(98) 



we have the so-called YZ piece, 



F'MU^') = ^KW, Ml Ml) + \n'i^^{q\ Ml Ml) 



6 



Ml Ml) + ^/Cl,Yz(g^ Ml Ml) - RsMq') 



(99) 



The quantity Rs^yz is rather complicated, so before writing it down explicitly we first develop 
some useful notation. For a quantity /(q'^, . . .), we define the auxiliary functions 



f{q\...)^f{q\...)-{q'-M')f{M\ 



(100) 



where the 'M^' quantities become Ml for the case of isospin fiavour and M^ for hypercharge 
fiavour. Then we have for i?3,YZ the expression 



R^,Yz{q^ 



q^-Ml 



p,M^-,MlMlMl) 



1 



) + -h,Yz{q';MlMl) 



+-/i,Yz(g'; Ml Ml- 2{Ml + Ml)) 

+ ^t/Yz(g'; M,^; Ml) + ^Ryziq'; Ml Ml Ml)) 
+l^Yz(g'; Ml Ml; 2{Ml + Ml)) 



1 . , . „ , 2(Ml - Ml) 
+ -RYz{q'; Ml Ml; ^ "3 



For example, we obtain at = the numerical value 

F2n^^^2(0) = 1.927 X 10"^ GeV^ 



;ioi) 



(102) 



B. The Spin-zero Isospin Polarization Amplitude 

For the isospin polarization function flf^ , we find for the remnant contribution. 



FltlfU^l') = ^ 



M^ 

71^ 



361 



+ 



49 C 



294912 36864 



J_ Ml 
3072 



11 



12288 



log 



Ml 



1 

— — log— ^ 
9216 ^ tJ? 



(103) 
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whereas the counterterm piece is given by 

For the piece coming from the finite functions, we have 

Flt\SU^') = 4^2,Yz(g^ Ml Ml) + \s,M<l\ Ml Ml 
+ ^52,Yz(g', Ml Ml) + i/C2,Yz(g', Ml Ml) 

wiinw.MiMi + Inw.MiM'^, 



6 



+l^nW. Ml Ml) + ^/Cl,yz(g^ Ml Ml) - R.^ql 



-2 



^/i,YZ {Ml Ml Ml Ml) + l/2,Yz {Ml Ml M 



+ ^/i,YZ {Ml Ml Ml; 2{Ml + M 

+ ^/i,YZ (m^^; M,^; Ml- \{Ml - Ml)) 
M'^ 

-^^Yz(M.'; Ml Ml) - -^SUMl Ml M^2) 
-^^Yz(M.^; Ml Ml Ml - ^R^Ml Ml M^ 2{Ml + Ml)) 
-^^^YziMl Ml Ml, \{Ml - Ml)) - \ul,^{Ml Ml Ml) 
The numerical value of = is found to be 



2rr(0) 



3,YZ 



(0) 



-2.954 X 10"^ GeV^ 



(104) 



(105) 



(106) 



(107) 



As mentioned earlier, the remaining finite results (hypercharge polarization amplitudes, 
etc) appear in Appendix B. 



C. Isospin Polarization Functions and MS Renormalization 

The discussion in Sect. VI-C allows one to re-express the isospin polarization functions 
il^l^^ in the MS renormalization. The result of this is, in essence, to replace the finite 
renormalization constants obtained in A-subtraction by the corresponding MS quantities 
and at the same time to omit all terms from the remnant contributions containing the 
constant C. 

As an example, let us determine the relation between the renormalization constants Pa 
and P^^. Starting from Eq. (p^ ) and applying the relations in Eq. (|90|), we find 

^'F = Pa + ^ (-2Bg> + bS') - g (-2Bg' + sf) . (108) 
Then from Eqs. (^) , (^3]) , we obtain 
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p, 



MS 



Pa- 



3 C 
16 (4^ 



This result is entirely consistent with the form obtained earlier in this section for fl 
is, from Eqs. (|95D,(]97|) we have 



(1) 

3 • 



F2 



-Pa + 



3C 



+ ... 



40967r4 

But this is just the combination of factors appearing in Eq. (|109|) and we conclude 



(109) 



That 



(110) 



F2 



pMS I 



;iii) 



Analogous steps lead to the further relations 

5 C 



VA 



dMS 



B 



MS 
ll,r 



pMS 

-°13,r 



32 (47r)4 
17 C 
96 (47r)4 

49 C 
' 576 (47r)4F2 
173 C 
5184 (47r)4F2 



;ii2) 



Again, the forms for flg^''''' obtained by us in A-subtraction are found to convert to MS 
renormalization by simply removing all C dependence from the polarization functions and 
employing the MS finite counterterms. 

Finally, we point out that for the renormalization constants P, Q, and R (c/Eq. (p3|)) 
which appeared in our two-loop analysis of vector-current propagators ^ , there is no differ- 
ence between the A-subtraction and MS schemes. This can be traced to the fact that the 
quantities P^^\ Q'"^\ and R^'^^ have only contributions from Lg^^ and L^q. 



VII. SPECTRAL FUNCTIONS 

As noted in Sect. II, there will generally exist two spectral functions, Pxai.'f') ^^^1 PAa(?^) 
for the system of axialvector propagators. In the following, we determine the 37r contribution 
to the spectral functions for the isospin case a = 3. The KKtt and KKrj components will 
have thresholds at higher energies. 



A. Three-pion Contribution to Isospin Spectral Functions 



Spectral functions can be determined from the imaginary parts of polarization functions 
;c/Eqs. (!),(!)), 

&{q') = -^rnli^\q') and pfa{q') = -^rnli^\q') . (113) 
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In two-loop ChPT, the imaginary parts arise solely from sunset graphs, and in terms of the 
notation established in App. A we find the Sir component of p^] to be 



1 



71 F3 fl67r2)2 



Xm 



(1) 



+ 



M: 



- - 4Z« 



(114) 



where the above finite functions are evaluated with 3tt mass values. There is a comparable, 
but rather more complicated, expression for p^g [Stt] which we do not display here. 



B. Unitarity Determination of p^3^^ [Svr] 

Unitarity provides an alternative determination of the three-pion component of the 
isospin spectral functions PA3[37r] and PA3[37r]. The first step is to relate the spectral func- 
tions to the fourier transform of a non-time-ordered product, 

piliq') W - q'gn + pflil') q'q' = ^ [ d'x e^'-(0|^^(:r)A^(0)|0) . (115) 

/vr J 

One obtains the three-pion contribution by simply inserting the Stt" and tt+tt^tt" interme- 
diate states in the above integral. To determine the relevant iS-matrix element, we employ 
the lowest order chiral langrian of Eq. (^ to find 

{as{q,\)\S\n'ip^)n\p,)7r%p,)) = 

{as{q,X)\S\7r^iPi)^-iP2)7i\p3)) = 

^(27r)^5W(g -PI-P2- P3)e;iq, X)M'^_o , 



;ii6) 



where the invariant amplitudes are given by 



Ml 



000 



v^F. ^ q' 



-0 



2p^o+q' 



Ml ' 
Ml 



2g -po 



Ml 



(117) 



Observe that A^+_o has two distinct contributions, a direct coupling and a pion pole term, 
whereas A^ooo has only a pion pole term. In the chiral limit of massless pions, the above 
amplitudes are conserved {q^M.^,-,, = 0) as required by chiral symmetry. 

To determine the spectral functions from Eq. (|115|) , we take p = = 3 for PasIStt] and 
p = 1/ = for PAsfSTr]. Throughout, we work in the Lorentz frame where = (g°, 0). With 
the tacit understanding that we consider only the three-pion component and take > 9M^ 
in the following, this leads to 
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J(l), 



1 



1287r6F2 
d^pi d^p2 d^P3 
El E2 E3 



{pirS^'\q - Pi - P2 - Ps) 



and to 



1 



Mi 



2 (0)/ 2\ _ 

j(o) ^ /■ d^d^(pp3 n_ 

J El E2 E^ \6 



/(o) 



(118) 



(119) 



- 4p!i + 4(p0)2^ ^(4)^^ -P1-P2- P3) 



The former experiences only the 7r+7r~7r° contribution whereas the latter has contributions 
from both the 7r"'"7r~7r° and Sir^ intermediate states. 

Integration of the above integrals is staightforward and one obtains a form involving 
integration over the invariant squared-mass of the 7111^2 subsystem. For example, for the 
spectral function pi^](g2) of Eq. (|Tl|) we obtain 



1 



da 



where 

A(x, z) = + + — 2xy — 2xz — 2yz 
With the change of variable 



^A3/2(g2^a,M2) 



a = - [^^q^ - 
we obtain the final form, 



2Mi 



1 



g2-M, -2M2 cos 



7687r4 F2 

where x = M^^j^fcp and /'^^^(a;) is the dimensionless integral 

x)= - (1 - x)^ - 2x2 / 
L2 J JO 

^ -4x2 + (i(l - x)2 + x2 - 1 + [i(l - Xf - 2x2 

X 



J{1), 



sin (1 + cos I 



a/2 



cos ( 



3/2 



i(l - X)2 + 2x2 + [i(l - X)2 - 2x2 



COS ( 



1/2 



In like manner, we find 



Mt 



5127r4 F2 (g2-M2)2 



7 



/r(a;)-4/r(x) + 2/^"^(x 



r(o) 



(0), 



e(g2 - 9M2) , 



(120) 



;i21) 



;i22) 



(123) 



(124) 



(125) 
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where 



1 - X? - 2x^ 



3/2 



X 



-4x2+ {^(1-xf + x^ -1 + 



sin 6 (1 + cos( 



-il-x? -2x'^ 



a/2 



cos < 



21 1/2 



(l - a;2 - i(l - a;)2 - - xf - 2x^] coscp) 



i(l-a;)2 + 2a;2 + 



i(l -x)2 -2a;2 



COS( 



1/2 



(126) 



The spectral functions obtained in the unitarity approach described here agree precisely 
with those obtained from the imaginary parts of the sunset amplitudes. 



VIII. CHIRAL SUM RULES 

In previous sections, we have determined the isospin and hypercharge axialvector prop- 
agators to two-loop order in ChPT. Essential to the success of this program is the renor- 
malization procedure by which the results are rendered finite. As a consequence of the 
renormalization paradigm, however, the physical results contain a number of undetermined 
finite counterterms. In particular, the real parts of the polarization functions 11^^ and 11^^ 
(a = 3,8) contain two such constants from one-loop order {L^q and h[^^) and five inde- 
pendent combinations of the {B^'^} from two-loop order. Of these, the constants Hi'^ and 
are related to contact terms which are regularization dependent and thus physically 
unobservable. However, the remaining counterterms (called 'low energy constants' or LEG) 
must be extracted from data. In this section, we describe how some of the 0{p^) countert- 
erm coupling constants are obtainable from chiral sum rules, and as an example we study a 
specific case involving broken SU{3). 

A derivation of the chiral sum rules together with an application of one of them appears 
in Ref . [^| . We refer the reader to that article for a general orientation. For the purpose of 



writing dispersion relations it suffices to note that at low energies the polarization functions 
are already determined from the results obtained in previous sections, although some care 
must be taken with the kinematic poles at = Q in the individual functions U^^l and H^l- 
As regards high energy behaviour, the large-s limit of spectral functions relevant to our 
analysis can be read off from the work in Refs. |2^J2^. However, since we are calculating 
up to order p^, terms up to and including quadratic dependence in the light quark masses 
must be included. Thus, for example, the asymptotic expansion for the isospin axialvector 
spectral function summed over spin-one and spin-zero reads 



„(l+0)/ N _ J_ 

PA3 ^')-s^2 





" 2" 

m 




1 + 12— 




1 . 



+ 0{all/s')] . (127) 



Similar expressions hold for the other components, and we summarize their collective leading 
asymptotic behaviour by 

P^ll{s)-0{1), p^^lis) ^ 0{s-') , {pil - P^lhis) - 0{s-') . (128) 
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The real parts of the polarization functions show exactly the same asymptotic behaviour 
as the imaginary parts, i. e. expressions analogous to Eq. (|128|) hold. The Kallen-Lehmann 
spectral representation of two-point functions then implies the following dispersion relations 
for the axialvector polarization functions of a given flavour a = 3, 8, 

g^nS(g^) - hm {q'utiq')) = ^ H ds -^^^ , (129) 



<j2=o ^ ^ <j2=o dg^ ^ ^ Jo 



s — — te) 



(130) 



We work with g^n^|'^^''(g^) due to the presence of = kinematic poles. Moreover, the 
subtraction constants have been placed on the left hand side in Eqs. ( |129|) ,( |T30|) in order to 
equate only physically observable quantities. Dispersion relations involving SU (3)-breaking 
combinations have an improved asymptotic behaviour, such as 



r(0) 



(1) 



r(0)^ 



ds 



(1) 



— le 



(131) 



and 



g^(nS - nS)(g^) - hm (g^ (nS - nS) (g^)) = g^ ds 



s — q"^ — ie 



(132) 



Sum rules are obtained by evaluating arbitrary derivatives of such relations at g^ = 0. 
For sum rules inferred from Eqs. ( |129| ),( P^ ) it is preferable to express the left hand side in 

Aa 1 



terms of n^°'^'* 



1 

n\ 



d 



n?l(o) 







{n-l)\ 



_d_ 



n n-1 



ni^l(o) - ^ 



_d_ 

'dUf 



{n > 1) 

nS(o) 



-ds^ 







{n > 2) 



(133) 
(134) 



where PaI{s) is defined in Eq. (^. Finally, Eq. ( |132|) leads directly to the following sequence 
of sum rules explicitly involving broken SU{3), 



1 



_d_ 

dq^ 



(nS + n 



(0) jtW jt{0) 
A3 " ^^AS " ^^A9, 



(0) 



{p^AL+pii-pii-pim 

on+l 



.(0) 



.(1) 



,(0)^ 



(135) 



where n > 0. 

For this last sum rule, let us consider in some detail the case = 0. An equivalent form, 
better suited for phenomenological analysis, is given by 



n 



(1) 

A3 



n 



d /^(o) 



(0) - (n,3 



n 



(0) 



°° , (pa3 - Pas 
ds 







(136) 
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Evaluation of the left hand side (LHS) of this sum rule yields 



1 f, A/f 2 — A/f 2 
LHS = — ' ^ QAif^) + 0.001053 



1 



F2ri67r2)2 



/i2 



log 2 




Ml f 5 1 , M|.~ 



+M2log-^^-- + -log 
+1287r^LSS log ^ - M^^ log 



(137) 



Recall from the discussion at the beginning of Sect. VI that there are three distinct sources 
for the finite low energy terms: (i) 0{p^) CTs {B'f^}, (ii) the 'remnant' contributions, and 
(iii) the finite Y, Z integrals (c/Eq. (0)). It is a combination of the latter two which give 
rise to the numerical term (which is scale-independent and vanishes in the SU{3) limit of 
equal masses) in the first line. We have displayed all chiral logs explicitly, and Qa(/^) is the 
0{p^) counterterm defined earlier in Eq. (0). Since the full expression is scale- independent, 
this allows one to directly read off the variation of the contributing counterterm combination 
at renormalization scale /x. 

We can use the sum rule of Eq. ( |137| ) to numerically estimate Qx{n). One needs to 
evaluate the spectral integral on the right hand side (RHS) of Eq. ( |136D . For our purposes, 
it is sufficient to approximate the contribution of the isospin spectral function in terms of the 
Oi resonance taken in narrow width approximation, p^p^^^^{s) ~ gaiS{s — M^J. Employing 
resonance parameters as obtained from the fit in Ref. [E^, we obtain 



ds ^ 0.0189 . (138) 

s 

Although consistency with QCD dictates that we also include the large-s continuum 0], the 
leading-order contributions would cancel in Eq. ( |136| ) and the remaining mass corrections are 
small. As regards the hypercharge spectral function p^g, little is presently known. The lowest 
lying resonances which contribute are /i (1285) and fi (1510) but the couplings of these 
resonances to the axialvector current have not been determined. Since the corresponding 
sum rule for vector current spectral functions exhibits large cancellations between the 
contributions from p (770), uj (782) and $ (1020), we expect a similar cancellation to be at 
work in the axialvector sector. To obtain a rough estimate, we assume the two resonances 
/i (1285) and /i (1510) can be approximated by a single effective resonance with spectral 
function p^s^is) ~ Qas^i^ — M^^). Assuming further (yf^g ~ and M^g ~ 1.4 GeV we 
estimate the hypercharge contribution to the RHS of Eq. (|136|) as 0.012. Allowing for a 
50 % error in this estimate places the RHS in the range 0.001 < RHS < 0.013 and leads 
finally to 

0.000043 < QA(MaJ < 0.000130 , (139) 



28 



where the renormahzation scale = Ma^ has been adopted. This is clearly to be taken as just 
a rough estimate. Only experimental determination of the missing coupling constants can 
provide a more reliable estimate. In addition, a more thorough phenomenological analysis 
will involve use of the entire spectrum. However, this example serves to illustrate the general 
procedure. 

We have not touched on chiral sum rules involving both vector and axialvector spectral 
functions. The most prominent example of this type is the Das-Mathur-Okubo (DM0) sum 



rule pO[ which, in modern terminology, has been employed to determine the LEG L^^^ 
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In Ref. 251 we have shown how the DM0 sum rule must be modified to be valid to second 



order in the light quark masses. Recently, r-decay data has renewed interest in this sum rule 



from the experimental side |31]. We have begun a phenomenological study of the DM0 sum 
rule using the two- loop results of polarization functions obtained both here and in Ref. 0. 
Results will be reported elsewhere. 

Finally, there are also those sum rules involving no 0{p^) counterterm coupling constants, 
i.e. those obtained by taking appropriately many derivatives of the dispersion relations 
Eqs. ( |129|) -( |T3^ ). From experience with the corresponding inverse moment sum rules of 
vector current spectral functions |T^ , we expect these sum rules in general not to be verified. 
This is because the relevant physics (which involves the low-lying resonances) enters the 
relations only in higher order of the chiral expansion. A quantitative study of these sum 
rules is deferred to a forthcoming publication. 



IX. CONCLUSIONS 

Our analysis of axialvector current propagators in two-loop ChPT has led to a complete 
two-loop renormahzation of the pion and eta masses and decay constants as well as the real- 
valued parts of the isospin and hypercharge polarization functions. It has yielded predictions 
for axialvector spectral functions and has allowed the derivation of spectral function sum 
rules. 

Despite the complexity of many of the individual steps and results, the sum of tree, 
one-loop and two-loop contributions to the axialvector propagator yields a simple overall 
structure, 

(140) 

where flavor labelling is suppressed. Comparing this to the general decomposition of Eq. (^ 
yields 



n«(g^) = 2Lffi - AH?^ + II«(g^) - El±Ip^ 



* " ' q^rf - AP) ^ ' 

As noted earlier, there are kinematic poles at = in both the spin-one and spin-zero 
polarization functions, but the sum 11^^ -|- 11^'' is free of such singularities. 
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A large number of 0{p^) counterterms entered the axialvector calculation, and many 
constraints among them were obtained from the subtraction procedure. Thus given the total 
of 23 0{p^) counterterms which appeared by employing the basis of Ref. |]2^, each of the 
and A subtractions were found to yield 14 constraints. The analysis of vector propagators in 
Ref. [0 yielded another 3 conditions for each of the A^ and A subtractions. This total of 17 
conditions constraining the {B^'^} and {B^^^} counterterms is of course universal and can 
be used together with results of other two-loop studies. We can summarize the remaining 
nine finite 0{p^) counterterms as 



Polarization Amplitudes 
Decay Constants 
Masses 



Pa,Qa,Ra, b[°i ,B[f 

{Be} (£ = 1,...4) (142) 
{Be} (£ = 1,...9) , 



where Pa,Qa,Ra are defined in Eq. (|9|) and the {B^} in Eqs. (|BToD , (p20|) . We have 
made preliminary numerical estimates for Qa in this paper (c/Eq. ( |139| )) and for Pa in 
Ref. p5[|. The counterterm B^^ is related to a contact term and is regularization dependent, 



much the same as the constant h[''> appearing in Eq. ( |140|) . However, these terms always 
drop out when physical observable quantities are considered. The constant Ra is seen 
to contribute equally to all flavour components of the axialvector polarization functions. It 
cannot therefore be accessed by the chiral sum rules involving broken SU (3) considered in the 
previous section. However, by combining the results obtained here with the two-loop analysis 
of the vector current two-point functions, the combination R — Ra is seen to constitute a 
mass correction to the Das-Mathur-Okubo sum rule. The details of this analysis will be 
presented elsewhere. Finally, little is known about the nine constants {Be} {i = 1, . . . , 9) 
which determine (together with the calculated loop contributions) the corrections to 
masses and decay constants. We have not attempted to estimate these constants {e.g. by 
the resonance saturation hypothesis) since as far as the axialvector two-point function is 
concerned, their contribution is implicit. However, the explicit expressions given here can 
be used in further ChPT studies when expressing bare masses and decay constants in terms 
of fully renormalized physical quantities. 

In view of the length and difficulty of the calculation, it is reassuring that a broad range 
of independent checks was available to gauge the correctness of our results. We list the most 
important of them here: 

1. Because the calculation involved independent determinations of isospin and hyper- 
charge channels at each stage, the SU{3) limit of equal masses provided numerous 
tests among the set of isospin and hypercharge decay constants, masses and polariza- 
tion functions. As a by-product, it also revealed the presence of previously unnoticed 
identities among the sunset amplitudes. 

2. As shown in Sect. VII, it was possible to determine spectral functions directly from 
the two-loop analysis or equivalently from a unitarity approach which employed one- 
loop amplitudes as input. In this way, both the specific sunset integrals as well as the 
structural relations of Eq. (|141|) were able to be tested successfully. 



3. It turned out that although most of the divergent terms could be subtracted away 
with counterterms, there occurs no counterterm contribution in M^f^"^ . To avoid 
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disaster, there must thus be a cancellation between sunset and nonsunset numerical 
terms. Such a cancellation indeed occurs and constitutes a nontrivial check on our 
determination of the sunset contribution. 

4. Given that n^''(g^) can be shown to vanish in the limit of zero quark mass, it follows 
that our final result in Eq. ( |14CI| ) has the correct chiral limit. 

5. Lastly, we have explicitly verified in MS renormalization that the constant C is absent, 
as must be the case. 

Yet more on this subject remains to be done. This is especially true of the material 
composing Sect. VIII, where an application of S'f/(3)-breaking sum rules to determine finite 
0{j)^) counterterms was discussed and additional points were raised. Future work will be 
needed to carefully analyze the axialvector sum rules, particularly the role of existing data to 
provide as precise a determination of the counterterms as experimental uncertainties allow. 
We can, of course, combine the results of the present axialvector study with the vector results 
of Ref. 10 to study an even wider range of sum rules {e.g. as with the proposed determination 
of i?A discussed above). On an even more ambitious level, our experience with such relations 
makes us optimistic about the possibility of describing a possible framework for interpreting 
chiral sum rules to arbitrary order in the chiral expansion. Finally, it will be of interest 
to reconsider the phenomenological extraction of spectral functions such as PassP^i"] and to 
stimulate experimental efforts to extract spectral function information involving non-pionic 
particles such as kaons and etas. 
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APPENDIX A: SUNSET INTEGRALS 



In the following, we compile mathematical details related to the sunset amplitudes which 
appear in the two- loop analysis. First, we give integral expressions for the quantities occur- 
ring in the IPI, vertex and self-energy amplitudes. Then we write down the finite sunset 
amplitudes which remain after the singular parts have been identified. Finally we discuss 
certain identities which relate various sunset integrals. Additional work on sunset integrals 
can be found in Refs. p2|JT0| , p!2|JT7[] for the equal mass case and in Refs. [p^- |35| , p!T| for the 
general mass case. 



1. Definitions of Sunset Integrals 



For the sunset amplitudes containing unequal masses, we shall denote the mass occurring 
twice as 'M' and the third mass as 'm' {e.g. for KKn amplitudes, we have M itlk and 
m — >• m^). The quantity Ti^i, appearing in the IPI sunset amplitude of Eq. p^ ) is defined 
by the integral expression 



Fin^,{q\m\M') = / dt 



f-2 _ ^2 



db 



1 



(62 - M2) • {{Q - 6)2 - M2) 



where Q = q — t and 



{S] S^] S^2j} — dt 



From covariance, we have 



,7 {^i^fii tfitu} 



^2 — 777.2 



db 



(62 - M2) . ((Q - 6)2 - M2) 



S,iq',m',M')^q,Siq',m',M') , 

S,,{q\m\M') = q^q,S,{q\m',M^) + g^,S2{q^m^,M^) , 



Also appearing in Eq. (|32|) is C^^, defined by 

1 



F2£^,(g2, 777^,M' 



dt 



j^2 _ jy^2 



db 



{Q-2b),{Q-2b), 
{b^ - M2) ■ ((g - 6)2 - M2) 



which can be expressed in the equivalent form 

1 



F^C^,{q',m\M 



d-1 

,2 _2 ]\4'2 



+ g^,ulC2{q\m\M')+Ag^ 



q^qylCiiq^^m^^M"^) 

'3 A-d" 



A{m^)A{M^ 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



where /Ci and /C2 are given respectively in Eqs. (|A18| ), (|A19| ) below. 

The vector-valued integrals Ji^ and which contribute to the axialvector vertex func- 
tion in Eq. (B3), are defined as 
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h, (f- m\ M\ a) ^ q^hW\ m\ M^- A) 



= / dt 



7 {q - 3t); 



■ ^^Q^ -2q-t + 2b-{Q-b)+A 



f^-m? J (62 - M2) . ((Q - 6)2 - M2) ' 



= / cit 



db 



{Q-2b),{Q-2b), 



(A6) 



(A7) 



t2 - m2 7 (62 - M2) ■ ((g - 6)2 - M2) ■ 

Finally, in the calculation of self-energies, there appear quantities S, R and U. The 
function S is already defined in Eq. (|A^), and we have for R and ?7, 



R{q^;m^;M^;A) : 
di 



t2 — 7722 

U{q'^;m'^;M^) = 
r di 



~ [jq -tf ~2q-t + 2b-{Q-h)+ K\' 
(62 - M2) ■ ((g - 6)2 - M2) 



^2 — m2 



c/6 



{{q + t)-{2b + t-q)y 



(A8) 



(A9) 



(62 - M2) ■ ((g - 6)2 - M2) ■ 

Analysis reveals that all the sunset contributions can be expressed in terms of the func- 
tions S, S, Si, 5*2, /Ci and IC2. This is already evident from Eq. ( |Al| ) for Hfj_u and from 
Eq. ( [A5| ) for One can deduce the additional relations 



Ji^ [q^; m^; M^; A) = -2q^ [2A\M^) + A{m^)A{M^) 

+ (g^S(g^ m^; M^) - 35^(g2; m^; M^)) [2(g2 + ^2 _ ^ ^ 



6 kg.5"(g'; m^; M^) - 3g-5^,(g2; m^; 



2F2 



2(3 + rf-4) ox ./,.2x 

'-A{m^)A{M^) 



(AlO) 



(All) 



as well as 



R{q^;m^;M^;A) = 

A^Siq""; m'; M^) - 4A [ - (g' + - M2)5(g2; m^; M') 

+ 3q^S^{q^; m^; M^) - A^{M^) + A{M^)A{m') 

+ 36g^g^5^^(g2; m^; M^) + 24(M2 - - q^)qf,S''{q^; m^; 
+ 4(g2 + - M^Siq^; m^; M^) + A2(M2)(4m2 - 12g2) 
+ A{m^)A{M'^){m'^ - 6g^ - Gm^) 
U{q^;m^;M'^) = 



(A12) 



rf- 1 



+ (g^ + m^) 7: + 



3 d-4^ 



A{m^)A{M^) 



+ q^ (g'/Ci(g2; m"; M^) + /C2(g'; m^; 



(A13) 
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The functions S, S, Si, S2, /Ci and JC2 can in turn each be written as the sum of 
terms (which diverge in the d ^ A hmit) proportional to gamma functions plus finite-valued 
functions {Y^^'''\q^ , m\ M^)} and {Z^''\q\m\ M^)}. Thus, we have 



and 
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(A14) 
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(A18) 
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/C2(g^m^M2) 
4{d 
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'(1)_2 



(A19) 



For the sake of simplicity, we have omitted the arguments of the {^J"-*} and the {Z^")}. 
We have verified that Eq. ( |A14|) agrees in the equal mass limit with the explicit expression 
appearing in Ref. [|17| . 



2. The Finite Sunset Integrals 

Having identified the singular parts of the sunset functions by expanding these quantities 
in a Laurent series about d = A, one can express the finite- valued functions {Y^f^} and {Z^^} 
which remain by means of integral representations, 

y ^-^ - TT.^ / , - 1 - ^ dxx- HI + Ag) , (A20) 

and 

r2\ 1/2+c 

r dn 

(16 



^i") = — ^ / da [1 dxx^ (ln(l + A^) - A^) (A21 

TT^)^ J4M2 \ a J Jo 



where 



Ag=[ g2 . (A22) 

\ X a 



For convenience, we shall introduce the dimensionless variables 

,2- l' „2_ ™' 



and likewise work with the reduced functions Yj"^^ and 



n) 



c 1 



Y^-\q\r^)^-^Y^-\q\r' 



167r2)2 ^ 
16^ 



Z^-\q'y)^r^,Z^:\<fy) (A24) 
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One is allowed to express such finite quantities in terms of the physical meson masses, and 
it is unerstood we do so in the remainder of this section. For the six flavour configurations 
which can contribute to the sunset amplitude, the parameter takes on the numerical 
values 



0.016 (ririn) 

0.020 {KKn) 

0.25 (Svr, 37]) 

0.31 {KKt]) 

3.82 (vrvrr/) . 



(A25) 



a. Behaviour at = and Near q"^ = 
In the = limit (i.e. m? = 0), analytic expressions can be obtained for Yj^^ and Z^"^ 
YtHf,0) = - t B(k+l;n+imk;a + V2) ^^^^^ 

k=l ^ 

and 

^,.,(,._ „) _ £ B(t+l; n+l)i?«=-l;a + 3/2) _ ^^^^^ 

k=2 ^ 

where B{m] n) denotes the Euler beta function. Observe in the summations that the indices 
begin at = 1 for FJ") and at A; = 2 for Z'^\ i.e. that 

FJ")(0,0) = Z(")(0,0) = Z(")'(0,0) = . (A28) 

For the more general case of nonzero but small g^, it is useful to employ a power series 

Y^-\fy) = FW(0,r2) +FW'(0,r2)g2 + ^Fj")" (0, r^)^^ + . . . 

Zi^\f,r^) = Zjr\0,r^) + Z(")'(0,r2)g2 + lz(")"(0, r^)^^ + . . . . (A29) 

For nonzero r^, one can obtain numerical values for the above = derivatives of 
r^) and Zj^"-\q'^,r'^). Of course, the integral representations of Eqs. (|A20|) , (|A21|) 
allow also for a straightforward numerical determination of the real part of the sunset am- 
plitudes for arbitrary g^. However, some care must be taken to obtain accurate values for 
g^ close to or above three-particle thresholds. 



b. Imaginary Parts 

For g^ < 1, the finite sunset amplitudes are real- valued. However, Vj^^ and have a 
branch point singularity at g^ = (1 + r)^ (corresponding to g^ = (2M + m)^) and become 
complex- valued for g^ > (1+r)^. We shall be concerned here with determining the imaginary 
parts of these quantities. 
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Consider first the integral X'^^\q^, r^) defined by 

XW(g2^r2)= j^dxx'' ln(l + A^) , 
wliicli can be rewritten as 



(A30) 



Jo 



dx 



dx x"- 



In ( a;' + a;( — - 1 - - ) + - 1 + ln{uf/x 



In {{x — x+){x — + Inluq^/x) 



(A31) 



where x± are given by 



^ ,^2 \ 2 /j^,^2 



(A32) 



The imaginary part of X'^")(g2,r2) will occur when the argument of the first logarithm 
in the above becomes negative, 

Im X^"'\<f jt"^) = / dx x^ Im ln((x — x+)(x — X-)) 
Jo 



Jx- 



dx 



so that 



Im Yt\q\ r') = /' -(1 - ( 

n + 1 ^ 



X 



n+l ry.n+1 



+ 



X 



(A33) 



(A34) 



where 



Uq 



(A35) 



The lower limit Uq on the li-integral is simply a reflection of the branch point occurring in 
the sunset amplitude at — (2M + m)^. Proceeding in like manner leads to the following 
formula for Xm Z^^^ , 



Tl ~\~ ± J Uq U ^ 



(A36) 
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3. Identities 



Given the set of sunset integrals {S; S^; S^,^}, it is not difficult to infer the following 
'trace identity', 

m\ M^) = m^S{q\ m\ M^) + A^{M^) , (A37) 

which is valid for arbitrary kinematics. 

It turns out that several more identities become derivable in the equal mass limit of 
S'f/(3) symmetry. This is a consequence of the symmetry constraint that the isospin and 
hypercharge results agree. Indeed, their direct comparison serves to check the correctness of 
the calculation. Interestingly, the identities discovered in the SU (3) limit are typically not 
at all a priori obvious. Below, we list and indicate the source of relations: 

1. Relating S to S: 

S{q\m\m^) = -S{q\m^,m^) . (A38) 
3 

2. IPI amplitudes: 

n^uiq'', m", m") = ^C^,{q\ m^, m^) . (A39) 



3. Vertex functions: 

hf,{q^] m^; m^; A) = m^; m^; 0) 

h^{q^- w?- w?- 0) = 3J2^(g2; m"- m^; 0) (A40) 



4. Self-energies: 



m^; m^; A) = R{q^] m^; m^; 0) + K^S{q^] m^; m^) 
i?(g2;m2;m2;0) = 3 f/(g2;m2;m2) . (A41) 



APPENDIX B: COMPENDIUM OF FINITE RESULTS 

In this Appendix we complete the compilation begun in Sect. VI of finite results in our 
two-loop calculation. We list in turn expressions for the hypercharge polarization functions, 
then the pion and eta decay constants and finally the pion and eta masses. The Appendix 
concludes with a brief summary relating our A-subtraction renormalization with the MS 



scheme of Ref. |jT7| 
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1. Hypercharge Polarization Functions 



First we display the corresponding hypercharge results, beginning with the remnant piece 



of the polarization function Ilg , 



F2fr(l) 

^vr S.rem 



Ml 



Ml 



77^ 



+ 



4096 
35 



r(l) 

13 
6144 

9 C 
1024 
3 C 



C I , Ml 

\ log — — 

512 512 ^ fi^ 



+ log 



Ml 

^2 



9 



37r2 



1024 



(0) 
^10 



1024 



log 



Mi; 



32768 4096 4096 
and then the counterterm contribution. 



3 Ml 
log^ 



F2 



F2 



R 



A 



F3 



Ra + -Qa 



(Bl) 



(B2) 



As was explained in Sect. VI, the l[q^^ dependence from the polarizations Ha^s.rem been 
removed via the definitions of Pa, Qa, -Ra given in Eq. (0). 
Concluding with the part from the finite functions, we have 



;,YZ 



\n\\{q\ Ml Ml) + \n\\{q\ MJ, MI) - Rs,yz 



where, making use of the notation established in Eq. (|100|) , 



R?,,Yz{q^ 



1 



+/i,Yz(g';M2;Mi 

(16M| - 7M2^2 



2(g2 - M2) 

2(3M|: - Ml 



h,Yz{q';MlMl: 



mi 



Ml) 



(g2 - M2)2 



M: 



YZl 



486 



^Yz(g^; M,2; Ml) + ip(g2; M^^ M^; ^ 



^ Mi; Mi 



+ \R{q'-Ml-Ml-'^-^ 



Ml 



and at g2 = 0, we find 



The last of the polarization functions is Ilg , for which the remnant piece is 



5.494 X 10"^ GeV^ 

(0) 



2971 



^ TT ■'■-'■8, rem 



Ml 



+ 



307 C 



497664 62208 256 



1 Ml 
log 



972 ^ /i2 



+ 



M2M2 / 8347 787 C 



+ 



25 



Mi 



TT^ V 995328 124416 4608 
/ 28123 2707 C 



log 



Ml 



7 Mf 
+ 7776 ^°^^, 



7962624 
9 



995328 



1 Mi 

og — 

3072 ^ /i2 



49 



M2 

log — — 

248832 ^ /i2 



4096 



log 



M|^ 



(B3) 



(B4) 
(B5) 



(B6) 
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Next comes the counterterm contribution, 



n 



(0) 



32 



?(o) 



8,CT — ^"tt \^ ~3~'^13 



'13 



32 



followed finally by the piece from the finite functions, 



9 



9 



F^fifU^') = lS,M^\ Ml Mi) + ^5,,Yz(g^ M|, Ml) 
+9'(^n^z(g^ Mi) + \n\W. Ml Ml) - i?s,Yz(g^)) 



-2 



^/i,YZ (m^2; M^; Mi; \{MI - M^) 



+\h,YZ [mi Ml Ml- ^(3Mi - Ml)) 



972 



12 

-^i^YzK; M^; Mi; ^(M,^ - Mi)) 
-^^YzK; Ml Ml- ^(3Mi - M^)) 

The numerical value of the YZ-part at = is 

F2fl^°^2(0) = -1.753 X 10-6 GeV^ . 



S[,^{MlMlMl 



(B7) 



(B8) 



(B9) 



2. Meson Decay Constants 

Before displaying explicit forms for -F^'^^], we must (i) implement the procedure (c/ 

Sect. VI) for removing all contributions from the {l\ counterterms and (ii) re-express all 
one-loop mass and decay constant corrections in terms of one-loop rcnormalized quantities. 

There are a priori seyen of the 0{p'^) counterterms l\ {£ = 1, . . . ,6,8) which contribute 
to the decay constants fI^^\-^ rem- Analogous to the procedure followed for the polarization 
functions n|3^8) rem' ^® rcmovc the {Lg~^^} dependence by defining the following four effective 
O(p^) counterterms, 

B, ^ (A<») - 3£;(»)) + f iS-" + f L<-' + f L<-'> 
B, ^ Fi - 2£;<°)) + |l<-'> + ?4-" + 54-" 
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5, 



26 



+ 



(-1) 



9 ' 9 



(-1) 
3 



(BIO) 



64-^^ 



9 

9 ^ 



9 

4-^44-^24-^^ 



It is natural to express the one-loop mass and decay constant corrections in terms of one- 
loop renormalized quantities. Thus, in the one-loop expressions for the decay constants and 
masses, we shall replace the tree-level parameters by their one-loop corrected counterparts. 



TT,K,7]) . 



(Bll) 



The quantities Ai,6i {i = n,K,ri) are compiled in Ref. (also see Eq. (pSj) of this paper). 
Thus, we write the decay constant relations for and through one-loop as 



Fo + — 



Ml , Ml 
^In^ 



167r2 



ALf (Ml + 2M|) + 4Lf 
+ . . . 



327r2 /i2 



ALf [Ml + 2M|) + 44°)m2 



3M 



327r2 



In 



M\ 



K 



+ ... 



(B12) 



and the mass relations for Ml, M]^ and through one- loop as 

Ml 



Ml 



Ml 



Ml 



+ ,r:-^ln 



-8 (4°^ -24°)) (m^ + 2M' 
M2 



327r2 



m\ + 
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- 8 (Lf - 24°0 Ml 
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- 8 (4°^ - 2Lr) {Ml + 2Ml) - 8 {if - 2L 



+ 



Ml Ml 
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M2 M2 
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2^ (3L(°) + L, 



(B13) 
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The expressions in Eqs. ( [B12| )- (|BT^ ) constitute our conventions for these quantities. In 
adopting our convention, we have made two kinds of choices: 

1. The prefactors I/-F0 (for decay constants) and 1/Fq (for masses) have been replaced 
respectively by and l/F^ since is the most accurately determined decay 
constant. 

2. We have employed explicitly throughout rather than use the GMO relation to 
replace it. 



A consequence of the procedure just described is to introduce additional 'spill-over' 
corrections at two-loop order to the set of decay constants and masses which were calculated 
earlier in the paper. Such contributions depend on the choice of convention discussed above. 
It is understood that the two-loop decay constants and masses (-Fji^rem; ^r)%m^ ^^Jemi ^Tvfrem) 
listed in the remainder of this appendix contain these spill-over corrections. Finally, for 
convenience we use the GMO relation in two-loop contributions to eliminate all factors of 
the eta mass not occurring inside logarithms. Any error thereby made would appear in 
higher orders. 

We have then for the pion decay constant. 
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and the corresponding counterterm amplitude is 

^'^i^CT = (2^1 - 2^2 - 4S4) - AMIMIB2 - SM^B^ . 
Finally, we find for the YZ contribution 

^'<YZ = ^/i,Yz(M,2; M,2; M,^; M,^; ) + 3^/i,yz(M^; M,^; M^; 2(M^ + M|,)) 

+1^ Vz(M,^ M^^ Ml- \{Ml - Ml)) + i/2,Yz(M^^ M^^; M^) 



(B15) 
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^S'^^{Ml Ml Ml) + ^5^z(M.^ M^^) 
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+\r!yz{MI; Ml; M^ Ml) + ^R'y^{Ml; M^ Ml; 2{Ml + Ml)) 
+^R'yz{Ml;M'^;Ml;^{Ml-Ml)) + ^U^^{MI; M^ Ml) 
The eta decay constant is treated analogously and one finds for the remnant contribution, 



i.3^(4) =M^(^^ 

n ^,rem n 1^ 1592524877^ 199065671^ 



5Lf - 84°^ - 8L(°) 



3 ^ 



+ 



14Lf 

-(0) 



64L?) 



24L 



(0)' 



^lio) _ ^ J- (0) 



, 72 ^ 144^^ 72' 



M 



M2 



+ -log^ 



49 5 , M| 

(0) 



8192 4096 
99 5 3 287r2 12 ^ 12 ^ 



(0) 



24^ 



(0) 



1 



log 



Ml 



24 ^ 40967r2^"'' /x2 



1 , M2 r 

+ — log ^ 



25 



+ 3^(0) + 3^0) 3^(0) 

122887r2 2 ^ 8 ^ 8 ^ 



9 



3443 



+MlMl 



(0) 



9 , Mi 



log 



40967r2 ° /x2 

127 C 



2488327r4 49766471^ 



44 



+ 324°) 



+ 



5L 



(0) 



8L, 



(0) 



(0) 



(0) 



(0) 



9 



18 



+ 



log 



Ml 



1 , Mi 
+ — log 



2 ^ 3 ^ +512^^°^^ 



41 



1 



- r(0) _ - r 



1 M? ^ 
+ — log 



,2 



184327r2 8 

187 



5^(0) 
24 



log 



20 487r2 ° PL 



2^(0) _ 2^(0) _ 1^(0) 

1244167r2 3 ^ 3 ^ 3 ^ 



1 



+M^^ ^^^^^ 



7303 C 



+ 



^Vl9906567r4 4976647r4 tt^ 



2r(0) 17 (0) 19 (0) 



+ 32Lf 



64 



+ ^4°) 7Lr + 2L, 



5Lf 

(0) 



8L, 



(0) 



-(0) 



1 , Mi 
+ — log^ 



11 



8L, 



+ 2Lf ) + -Lo 



(0) 



8L 



(0) 



8L 



(0) 



(0) 



-L 



(0) 



5r(0) 

4^4 



1 



_ 2r(0) , 

3^' +5127r2 



log 



Mi 



+ — log ^ 



211 



62 2 087r2 3 ^ 



3^3^ 



(0) 



-L 



(0) 



1 



log 



3 " 2567r2 - /X- 

whereas the counterterm and YZ amplitudes are respectively 
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3. Meson Masses 



In order to remove the seven l\ 0{p'^) counterterms 
define the following five effective 0{p^) counterterms, 
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Beginning with the pion squarcd-mass, we have 
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-^i?Yz(M^; Ml; Ml; '^-{Ml - Ml)) 
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4. Relation to MS Renormalization 

Our formulae for the hypercharge polarization functions flg^'^'' displayed earlier in this 
appendix refer to the A renormalization used throughout this paper. To obtain corresponding 



expressions in the MS renormalization of Ref. IJl^, it suffices to follow the discussion in 
Sect. VI-C for the isospin polarization functions. 

However some additional analysis is required in order to compare our A-subtracted decay 
constant and mass formulae to the MS scheme. We shall omit detailed derivation and simply 



50 



display the results, as the procedure mirrors that used to obtain Eq. (jll2|) . What is needed 
is the following set of relations between the A-subtracted constants {Bi} and the associated 
WS quantities Bf^ for £ = 1, . . . , 9, 



Be = Bf^-ABe , (B27) 



where 



AB (-— ^ ^ ^ ^ -1^ ^ -1] (B28) 

^ (47r)4 V 576' 96' 5184' 192' 27' 81' 288' 1296' 32/ ^ ^ 

As expected from our previous discussion, to obtain the masses and decay constants in the 
MS renormalization of Ref. one needs only make the replacements Bi —>■ Bf^^ and omit 



all dependence on the constant C. 



Figure Captions 



Fig 


. 1 


Lowest-order graphs for the axialvector propagator. 


Fig 


. 2 


One-loop graphs. 


Fig 


. 3 


Generic corrections to the axialvector propagator. 


Fig 


. 4 


Two-loop IPI non-sunset graphs. 


Fig 


. 5 


The two-loop IPI sunset graph. 


Fig 


. 6 


Two-loop IPI vertex graphs. 


Fig 


. 7 


Two-loop IPI self-energy graphs. 
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